Static versus dynamic friction: The role of coherence 
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A simple model for solid friction is analyzed. It is based on tangential springs 
representing interlocked asperities of the surfaces in contact. Each spring is given 
a maximal strain according to a probability distribution. At their maximal strain 
the springs break irreversibly. Initially all springs are assumed to have zero strain, 
because at static contact local elastic stresses are expected to relax. Relative tan- 
gential motion of the two solids leads to a loss of coherence of the initial state: The 
springs get out of phase due to differences in their sizes. This mechanism alone is 
shown to lead to a difference between static and dynamic friction forces already. We 
find that in this case the ratio of the static and dynamic coefficients decreases with 
increasing relative width of the probability distribution, and has a lower bound of 1 
and an upper bound of 2. 

While the facts that dry solid friction is proportional to the normal load at the contact 
and does not depend on the apparent contact area were established experimentally at least 
as early as in the 16th century by Leonardo da Vinci and are now known under the names of 



Amonton (1699) or Coulomb (1781) 
between static and dynamic friction 



it was probably Euler (1750) who first distinguished 
This difference has been explained in several, con- 



ceptually different ways. The reason was identified as: A collective depinning phenomenon 
0, the time strengthening of individual pinning sites 0,0], the shear melting of a lubrication 
ilm fj], mobile impurities at the interface |7], or the formation and healing of microcracks 
^]. The fact that all these mechanisms lead to the same macroscopic phenomenology raises 
the question whether they can be classified in terms of more abstract concepts. 

An attempt in this direction was made by Caroli and Nozieres [9j|, who proposed a model 
for dry solid friction based on the following physical picture: The surfaces have randomly 
distributed asperities which get interlocked. These interlocked asperities act as pinning sites 
resisting tangential motion. Under tangential load they are deformed up to a threshold 
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FIG. 1: a) Schematic plot of the force F as a function of displacement s is presented. At the 
threshold i the spring breaks and immediately reattaches, b) By reversing the displacement 
the force also changes sign. During the complete cycle O — > A ^ B — > C — > O energy is 
dissipated irreversibly and is numerically equal to the shaded area. The size of the hysteresis loop 
is proportional to the threshold length I, 

(which these authors call "spinodal limit"), where they break irreversibly releasing their 
energy in the form of phonons into the bulk. The threshold is a measure for the pinning 
strength. They argue that their model does not lead to a difference between static and 
dynamic friction, unless the strain of the pinning sites has different statistics in the static 
and in the sliding case or aging is taken into account. The latter aspect has been further 
investigated in ref. \^ and explains also the experimentally observed time strengthening and 
velocity weakening of the pinning sites. 

However, time strengthening is a slow process. This motivates us to explore in more 
detail what would be the influence of strain statistics at the pinning sites on the static and 
dynamic friction coefficients, /i s and \x&. Although time strengthening will not be considered, 
it can be included in addition to account for, e.g., velocity weakening. 

The model we consider in the following captures, we believe, the essence of the physical 
picture described above and at the same time highlights the concept of coherence, which 
is an ingredient in several different models (see (lOj and references therein). For the sake 
of clarity and analytical tractability we work out only a one- dimensional version, but the 
extension to the two-dimensional case is straightforward. 

Consider two solid bodies in contact, one being the fixed substrate (the "track") and the 
other the body to be displaced (the "slider"). To keep the equations simple, we consider 
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only motion in a fixed direction. Reversing the direction would lead to hysteretic behaviour 
like in Fig. ^ The friction force arises from interlocked asperities in the contact area, which 
are modelled here by linear springs with zero equilibrium length (see also U|), which only 
act in the tangential direction. Each spring has one end attached to the slider, while the 
other end is attached to the track. When the slider moves, each spring gets stretched up to 
an individual threshold length £, where it breaks. The elastic energy stored in a spring is 
completely dissipated when it breaks. 

In contrast to previous work we take here explicitly into account that the interlocked 
asperities are characterized by different threshold lengths £ with a probability distribution 
p(£), normalized as p(£) dl = 1. There is experimental evidence that this distribution is 
approximately Gaussian centered around a characteristic length 12]. 

In general the number of pinning sites and the distribution of their strength p(£) will 
change with time during a transient until steady state sliding is reached. However, one can 
imagine experimental situations where both are time independent, at least in an average 
sense Here we make this assumption deliberately in order to show that a difference 

between static and dynamic friction can ensue, even if the number of pinning sites and the 
distribution of their strength are time independent. With this assumption a new spring 
with the same parameter £ has to become active whenever one breaks. Hence the elastic 
restoring force from springs of threshold length £ becomes simply a sawtooth-like function of 
the displacement s (see Fig^). This displacement is assumed to be the same for all springs 
(approximation of a rigid slider). The friction force is the sum of all these elastic restoring 
forces. In the following all spring constants k are assumed to be the same, but this is not 
crucial. For example we checked that spring constants proportional to £ give qualitatively 
the same results. 

The crucial ingredient of our model will be discussed now. During sliding all springs 
will be stretched in the sliding direction by a random fraction of their threshold lengths. 
When the motion is stopped, the slider will recoil so that some springs get stretched in the 
opposite direction until the net force on the slider is zero (in the absence of an external 
shear force). In contrast to [9| we assume here that then the strain distribution becomes 
narrower, because the springs have time to relax. We have to discard plastic flow as the main 
relaxation mechanism when relative motion comes to a halt, since the speed of this process is 
proportional to the difference between applied stress and yield stress and should be therefore 
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too slow. A different mechanism is required, which is slow compared to the life time l/v of 
stretched springs during sliding with velocity v, but fast compared to the time a stationary 
contact is at rest. One possibility might be viscoelasticity: During sliding a nonequilibrium 
density of point defects in the immediate neighborhood of the surface is created. These 
point defects can be viewed as a viscous fluid penetrating the crystal lattice: if the lattice 
is exposed to time dependent stresses at very high frequencies during sliding, the defects 
hardly have time to diffuse and contribute to stress relaxation. However, if the frequency is 
low or even zero, as in the static case, the point defects move due to thermal activation to 
regions where they reduce the elastic energy of the entangled asperities (we note that this 
is different from plastic flow, which is due to the motion of dislocations). Diffusion of point 
defects is also a slow process, but since the distances are at the nm to /im scale, and since it 
may be assisted by strain, which can considerably reduce activation energies, it is still faster 
than plastic deformation. Consequently we expect that microscopic interface strains relax 
relatively fast in the static, but not fast enough in the sliding case. Conceptually this is 
different from time strengthening, where atomic diffusion would shift the threshold lengths 
£ towards larger values, an effect that occurs in addition, but is neglected here for the sake 
of working out the effect of stain coherence more clearly. 

In the context of our model we actually consider the extreme case, where all springs 
relax to zero elastic energy, as soon as sliding stops (full "coherence"). With the assumption 
that all springs are relaxed initially, the friction force as a function of displacement s for an 
apparent (macroscopic) contact area A and the number density n p of active pinning sites 
(or springs) reads 

oo 

f(s)=An p k Jp(£)t(£,s) d£, (1) 
o 

where t(£, s) = s mod £ is a sawtooth-shaped function of periodicity £. The phase of this 
periodic function is </>(£, s) = t(£, s)/£, which is a number between and 1. The behaviour 
of eq. ((H) is closely related to the probability distributions of these phases, 

w{4>, 8 ) = Jp(*)tf^-*M) dL ( 2 ) 

For example w(0, s) is the probability density of springs that break at displacement s. As 
these springs still contribute their elastic restoring force to f{s — ds/2) but no longer to 
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f(s + ds/2), the derivative of eq. (JTJ) is given by 

f'(s)=An p k [l-w(0,s)}. (3) 

This can be derived from eq. using the expression ((SJ) given below. 

Whereas the initial state is coherent in the sense that w((f>, 0) = 5 ((f)), coherence gets lost 
for large displacement, where all phases become equally likely: 

lim w(d>, s) = 1. (4) 

s— >oo v * 

To prove this we evaluate the integral in eq. (J2J) around each of the discrete values of £ for 
which the argument of the 5-function vanishes and obtain 

^s)=±p(^—)-^—. (5) 

^ \m + (j)J (m + (f)) 2 

Introducing the variable x = (m + (f))/s, which becomes quasi-continuous for large s, this 
converges to the Riemann-integral (provided p(£) is a Riemann-integrable function) 

oo ^ oo 

lim s) = J p (-) — = y = 1, (6) 



where the variable transformation y = 1/x was used, and the last equality is just the 
normalization of the distribution. This shows that after sufficient displacement the system 
forgets its initially coherent state. An important consequence of this decoherence is that 
the friction force for large displacements becomes constant. This follows immediately from 
eq. Q, because w tends to 1. Then the value of the friction force in eq. (0) coincides with 
its average, 

(f) = An p -jp(£)£ d£ = An p -(£), (7) 
o 

where {£) is the average maximum spring length. 

Figure shows the friction force f(s) as a function of displacement s on a microscopic 
length scale, for three different distributions p(£). For the numerical evaluation it was useful 
to rewrite the force f(s) in eq. ((H) by integrating eq. ©, after formula (jHJ) has been inserted: 

s/m 

I P(s) s ds 



f(s) = An p k 



* — 

m=l 







(8) 



In this expression the integrals can be calculated analytically for simple distributions, e.g., 
truncated Gaussians. All curves in fig. have in common that they start from zero with 
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FIG. 2: a) Friction force / as a function of displacement s of the slider for different maximum 
spring length distributions p(i). (/) is the dynamic friction force, £ c is a characteristic maximum 
spring length. Solid line for a truncated Gaussian distribution, p(£) cx exp[— {£ — £ c ) 2 /2(A£) 2 ] with 
width A£ = 0.15 t c . Dashed line for the same Gaussian distribution, but with width A£ = 0A£ C . 
Dashed-dotted line for p(£) cx 1/[1 + (£/Q% width is A£ = oo. Note that the height of the 
first peak decreases with increasing width, b) Friction force / as a function of displacement s 
of the slider for p{£) cx exp[— £ 2 /2(A£) 2 ]. The maximum value of / is approximately 1. 098(7} a t 
displacement s ~ 1.087A£. 

a slope of 2 in the natural units, (/) and £ c , and converge to 1 for large displacements. 
For the truncated Gaussian distributions the force has a maximum at a displacement close 
to £ c , where a large number of relatively strong springs are just about to break. As long 
as an external driving force remains smaller than this maximum, the position of the slider 
will move only on the scale of £ c . Therefore we interpret the first peak of the force as the 
static friction force, the pulling force needed to initiate motion of the slider. Once the object 
moves, the force necessary to maintain its motion decays to a smaller value with damped 
oscillations, which die out again on the scale of £ c . Therefore we interpret the asymptotic 
force as the dynamic friction force. 

We found numerically that the maximum is less pronounced the wider the Gaussian 
distribution is for a given £ c . This means that the difference between static and dynamic 
friction decreases. An interesting question is what happens for a distribution with finite 
average value, but infinite width. An example is shown in fig. |2K- in this case the force 
monotonically increases to (/), so that according to our interpretation the static and dynamic 
friction coefficients are equal. However, we regard this example purely as an illustration. As 
mentioned above, the empirically found distributions are approximately Gaussian. 
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At this point we can clearly formulate the main message of this paper: In our model the 
presence of an initial peak of the friction force, meaning that the maximum static friction is 
larger than the dynamic friction, is the result of an initial coherence in the strain distribution 
of the interlocked asperities. The height of the peak depends on the distribution of threshold 
lengths (or loop sizes, see FigHJ. After displacement of the order of few times the average 
threshold length the initial coherence is forgotten, the strains get out of phase, and as a 
consequence the friction force decays to the dynamical friction force. While we assume a 
complete initial coherence in our model, this is not absolutely necessary: Decreasing the level 
of initial coherence still results in a peak of the friction force, although with a decreased 
height (results not shown here). 

Now we show that the ratio of static and dynamical friction coefficients does not change 
under a rescaling of loop sizes £ — > a£. With the probability distribution 

p{£) -> ap(a£), (9) 

one gets according to eq. (JTJ) that the elastic restoring force at displacement s and hence 
also its maximum value (static friction force) and its average value (dynamical friction force) 
transform as 

f(s) -> -J(as). (10) 

Therefore the ratio between static and dynamic friction remains invariant. 

As mentioned above, the probability distribution p{£) is approximately Gaussian in many 
cases [3], i.e. completely characterized by its first and second moments, (£) and (£ 2 ). As 
under a rescaling £ —>■ a£ the first and second moments scale differently, {£) — > a(£) and 
(£ 2 ) — > a 2 (£ 2 ), but the ratio /x s //Xd remains invariant, it cannot depend on (£) and (£ 2 ) 
separately, but only on the invariant combination (£ 2 )/(£) 2 : 

The numerical analysis shows that g is a decreasing function of its argument (this tendency 
can be seen in fig. EtO- Therefore the extreme case, where the width of p(£) is zero \p{£) = 
5(£ — £ c )} gives an upper bound on the ratio of the friction coefficients. This case is special, 
as coherence never gets lost: All springs get stretched up to £ c , break simultaneously and 
are replaced by fresh, unstrained springs, which again get stretched up to £ c and so forth. 
(The mathematical reason why our proof that the force converges to the average value is 
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not valid in this case is that the S function is not Riemann-integrable.) Because of the lack 
of decoherence, the model gives stick-slip motion. However, for dry solid friction this is 
unphysical: Any small randomness for example of the times, when new interlockings form, 
would have the effect, that the springs ultimately get out of phase. Then the sawtooth- 
oscillations of the force f(s) = An p kt(£ c , s) would be damped similar to the oscillations 
shown in fig. E^l and would converge to the average value (/) = An p k£ c /2, which we still 
identify with the dynamic friction. Obviously it is half the maximal value of f(s), which 
we identify with the static friction. This gives an upper bound of 2 for the ratio between 
the friction coefficients in eq. (fTTJ) . Note that the upper bound would increase if time 
strengthening were taken into account. 

Together with the lower bound obtained if the relative width of p(£) tends towards infinity, 
we conclude that the model restricts the ratio of the friction coefficients to the interval 



Actually, for Gaussian distributions, a more stringent lower bound, approximately 1.098, can 
be given. This value is obtained, when the argument of g in eq. (|TT|) tends to infinity, which 
corresponds to the limit A£/£ c — > oo. Then the force f(s) approaches the one obtained for 
£ c = 0, i.e., for p(£) oc exp[— £ 2 /2(A£) 2 ]. For this case the ratio of the friction coefficients is 
approximately 1.098 (see fig. |2t>) independent of A£, due to eq. (JTUJ). Lower values can be 
obtained if p(£) has a power law tail, e.g., p{£) oc 1/[1 + (£/£ c ) 3 ], as shown in fig. EK- 

There is another interesting consequence of this theory. According to the theory of 
Bowden and Tabor the microscopic contact area of the pinning sites adjusts quickly by 
plastic flow such that the local stress drops to the yield threshold a c of the material. Then 
the normal load is equal to the real contact area times a c : 



where we assumed that a pinning site of loop size £ contributes a£ 2 with a constant geometry 
factor a of order 1 to the real contact area. Combining this with eq. (J7J) 



1 < ^ < 2. 



(12) 



/ n = An p cr c a(£ 2 ), 



(13) 



ft 



(14) 



one finds that 




(15) 
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As the friction coefficients should be independent of / n , we conclude that the number of 
pinning sites, An p , increases proportional to the normal load. This argument is not entirely 
compelling, as the pinning strength £ needs not be directly related to the microscopic contact 
area, and also the spring constant k might depend on the normal force. 

In this work we presented a simple model of dry friction, which explains why static friction 
force can be larger than dynamic friction force, in terms of the concept of coherence. 
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